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We present an improved version of the Stillinger-David polariza- 
tion potential of the intermolecular interaction in water. A clear 
algorithm of construction of a function describing the oxygen- 
hydrogen interaction in water molecules is formulated. A new 
approach to the modeling of a function screening the charge-dipole 
interaction on small distances is developed. To describe the long- 
range asymptotics of the intermolecular potential, the bare Stillin- 
ger— David potential is supplemented by a term related to the inter- 
action of dipole moments of oxygen ions. In addition, we introduce 
a term involving a deformation of the electron shells of oxygen 
ions to the polarization component. These corrections allow us to 
successfully reproduce all essential results of quantum mechanical 
calculations of the interaction energy for water molecules obtained 
by Clementi. Analyzing the behavior of the dipole moment of a 
water molecule as a function of the intermolecular distance, we 
obtain the estimate of irreducible two-particle effects in water. 



1. Introduction 

Astonishing properties of water are directly related 
to the formation of hydrogen bonds between water 
molecules [IHl]- Various potentials are used to simulate 
the influence of hydrogen bonds, as well as electrostatic 
multipole interactions between water molecules [IHH]. 
One has to distinguish between rather simplified poten- 
tials, which mainly describe the interaction between wa- 
ter molecules [T^TUl , and potentials, which describe both 
the molecule-to-molecule interaction and the interaction 
between model charges inside a water molecule |¥H51[T^ . 
As a result, the calculation, e.g., of vibrational spectra 
of water molecules becomes possible. 

In the framework of the Malenkov-Grokhlina-Poltev 
potential model [TT], a water molecule is simulated as 
that composed of three effective charges: one positive 
and two negative charges located at the centers of the 
oxygen and two hydrogen atoms, respectively. The Jor- 
gensen potential |9j also represents a water molecule 
as three effective charges: the positive charges are 
centered at hydrogen atoms, but the negative charge 
is displaced with respect to the oxygen anion center. 
It was done to match the dipole moment of a water 
molecule. The model of intermolecular interaction po- 



tential turned out optimal, provided that the charge 
values are fractional. Note that the Jorgensen and 
Malenkov-Grokhlina-Poltev potentials, owing to their 
simple structure, are widely used, while simulating the 
behavior of aqueous systems by methods of molecular 
dynamics. 

The Stillinger-David (SD) potential [ni[T3] is one of 
the most fruitful models of interaction between parti- 
cles, which was used to describe the formation of hydro- 
gen bonds in water. In the framework of this potential, 
a water molecule is represented by three charges, and 
the oxygen anion is supposed to be polarizable. With 
the help of the SD potential, both the intermolecular in- 
teraction and the vibrational spectra of water molecules 
are reproduced satisfactorily |H[T ^ IT ^ . The approach by 
Stillinger and David allows one to obtain a good agree- 
ment with experiment for the dipole moment of a water 
molecule and the angular dependence of the interaction 
between molecules which corresponds to the formation of 
hydrogen bonds. This potential was also used, while con- 
structing the autocorrelation functions for translational 
and angular velocities of water molecules, as well as when 
calculating their self-diffusion coefficient The SD po- 
tential was also used in works [rH - fTf] . where the spectra 
of acoustical and optical excitations in crystalline ice and 
the ice dielectric permittivity were calculated |17| . 

In works f¥'-'S], a modification of the Stillinger-David 
potential-the modified polarization model (MPM) - was 
proposed. This model substantially simplifies the form of 
potentials for the oxygen-hydrogen and oxygen-oxygen 
interactions and uses only one screening function, S{r), 
rather than two, 1 — L(r) and 1 — K{r), as the SD poten- 
tial does. As a result, the values of force constants were 
calculated more precisely, and the computation proce- 
dure was made simpler [3]. 

At the same time, the results obtained by Stillinger 
and David, as well as by the authors of MPM do not 
agree with the results of quantum-chemical calculations 
carried out by Clementi [TD] at distances shorter than 
about 2.5 A. A comparison of MPM and SD poten- 
tials with the interaction energy that was determined 
in work |10| evidences their unsatisfactory agreement at 



ISSN 2071-0186. Ukr. J. Phys. 2011. Vol. 56, No. 3 



225 



I.V. ZHYGANIUK 



distances that either correspond to or are much longer 
than those corresponding to the dimer formation. 

In this work, a generahzed version of the StiUinger- 
David potential (GSD) has been developed. In this 
model, we try to combine the most successful features of 
the SD and MPM potentials. In particular, (i) the func- 
tions that describe the energies of the oxygen-hydrogen 
and oxygen-oxygen interactions are similar to the ex- 
pressions obtained in the MPM; (ii) in order to deter- 
mine the function that describes the oxygen-hydrogen 
interaction energy in a water molecule, a very rigorous 
algorithm has been formulated. It takes into account the 
fact that the solutions of the system of algebraic equa- 
tions for the coefficients of a polynomial that is used to 
approximate the screening function 1 — L{r) have a very 
irregular character; (iii) besides the Coulomb and polar- 
ization components of the SD potential, the interaction 
between the induced dipole moments of oxygen atoms 
is taken into consideration, which is very important for 
a correct reproduction of the dipole-dipole interaction 
energy between two water molecules separated by long 
enough distances. As shown below, the screening func- 
tion 1 — L{r) changes monotonously, thus being substan- 
tially different from that given in works [l^.TS*. Owing 
to the aforesaid improvements, the results of quantum- 
chemical calculations by Clementi were reproduced quite 
successfully. 

In Section 2 of the work, a modified version of the 
Stillinger-David (MSD) potential is presented. It is used 
to obtain the dipole-dipole asymptotics of the interaction 
energy between water molecules. In Section 3, the laws 
of interaction between the hydrogen atoms and the oxy- 
gen one, which are included into a water molecule, are 
generalized. A short discussion of the results obtained is 
presented at the end of the article. 

2. Modified Stillinger-David Potential 

In the framework of the Stillinger-David potential ap- 
proximation [12J , the energy of interaction between two 
water molecules is expressed by the formula 



where the subscripts i,j = 0,1,2. The subscript i enu- 
merates the charges of one water molecule {i — corre- 
sponds to the oxygen charge, whereas i = 1, 2 to the hy- 
drogen ones), and j enumerates the charges of the other 
molecule. The charges are measured in terms of elemen- 
tary charge units, so that Qq — —2 and qi = (72 = 1- 
The Coulomb interaction between hydrogen and oxygen 
charges at large distances is reduced to the interaction 
between the dipole moments fi^ of water molecules as- 
sociated with the charges of hydrogen atoms in those 
molecules. 

The second contribution $n corresponds to the repul- 
sion potential of a hydrogen atom from the electron shell 
of the oxygen atom. It is approximated by the Born ex- 
ponential dependence 



$11 = 61 



E 



E 



-P2l'Oi02 



''O1O2 



(3) 



where bi and 62 are the amplitudes of hydrogen-oxygen 
and oxygen-oxygen, respectively, repulsion energies; pi is 
the reciprocal range of action of repulsion forces between 
the hydrogen atoms and the electron shell of oxygen; and 
P2 is the reciprocal range of action of the repulsion forces 
between the electron shells of oxygen atoms. 

The third contribution $in in Eq. (1) is responsible 
for the potential of interaction between the point-like 
charges of the first molecule with the polarizable oxygen 
atom of the second one. The oxygen atom is polarized 
under the influence of a field created by the charges of 
hydrogen atoms belonging to the same water molecule, 
as well as by the charges of hydrogen atoms and the po- 
larized oxygen atom belonging to the second molecule. 
The oxygen polarization gives rise to an emergence of the 
dipole moment /Xq, which characterizes the deformation 
of oxygen electron shells. Hence, the polarization con- 
tribution <i>in can be written down as follows: 



$ = $1 + $1 



III- 



(1) 



$111 = 



'■O1O2 



^ [1 ~ L{roi02 



0102 



The first contribution $1 in Eq. (1) is associated with 
a direct Coulomb interaction between oxygen and hy- 
drogen atoms in a water molecule, 



(2) 



'"01O2 
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(4) 



where //qi i and /Xqj are the dipole moments of oxygen 
atoms of the first and the second, respectively, water 
molecule; and 1 — L{r) is the screening function (see be- 
low). The dipole moments are determined in the molec- 
ular coordinate system (MCS), the origin of which co- 
incides with the center of mass of the oxygen atom in 
the water molecule. At large distances, the polarization 
contribution $ni is reduced to the interaction between 
the polarized oxygen atom of the first molecule and the 
charges of the second molecule. 

One can verify that, at distances that exceed the di- 
mensions of water molecules very much, the interaction 
potential $ looks like 

* = 'J'd(/^1H, M2h) + *d(MlH, Atos) + *d(M2H, /^Oi) , (5) 

where 

M2) 3(/xi 



*d(Ml,M2) 



01O2 



(M2 • roaOi) 



(6) 



In formula (5) , //jh and the contributions to the 

dipole moment of water molecules resulting from the spa- 
tial axrangement of hydrogen charges. The oxygens are 
at the origin of the MCS, and their contributions /Xo^ 
and /X02 to the dipole moment may arise only due to 
their polarization by the electric field of hydrogen atoms. 
Asymptotics (5) is incorrect, because the energy of in- 
teraction between water molecules at large enough dis- 
tances is determined by their total dipole moments. 



(7) 



We note that the infiuence of hydrogen atoms that 
belonging to the same molecule is by two orders of mag- 
nitude larger than the influence of charges in the second 
molecule. This fact completely agrees with the order of 
magnitude of a relative variation in the frequency of va- 
lence vibrations in the water molecule, when changing 
from the liquid- water to the saturated- vapor state LlSj . 

Unfortunately, asymptotics (5) is also characteristic 
for the modified polarization model (MPM) that was 
developed in works [4 -6J. 

To correct the asymptotics of formula (1), we have to 
pass to the modified Stillinger-David potential. 



$ = $7 + $11 + <f> 
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Fig. 1. (a) Configuration of molecules and (b) the angular depen- 
dence (at a = —360° -f^ 360°) of the interaction energy between 
water molecules calculated with the use of the SD (dotted curve) 
and MSD (dashed curve) potentials and the Clementi function 
(solid curve) 



where the component <I>iv looks like 
$rf(/Xoi,/Xo2)[l - K(ro^o^la)\. 



IV 



(9) 



Here, 1/a is a coefficient that modifies the screen- 
ing function 1 — K(r) for the interaction between the 
dipole moments of oxygen atoms in two neighbor water 
molecules. The optimal fitting of dementi's data (see 
Fig.[T]) was obtained at a = 2.235. 

It should be noted that the component <I>iv (Eq.(9)) 
of the MSD potential supplements the SD potential 
by including the dipole-dipole interaction between po- 
larized oxygens in water molecules. The effect of the 
dipole moment screening for oxygen atoms in water 
molecules is described by the function 1 — K(r ja). The 
screening function 1 — K(r ja) is similar to the function 
1 — K(r) which was used in works |12[[T3| to describe the 
dipolecharge interaction in the SD potential (see below) . 

Owing to the screening functions 1— i4r(r/a), the inter- 
action energy between water molecules almost does not 
change at short distances. Meanwhile, at large disances. 
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formula (8) gives the correct asymptotics, 

$ <I>d(/Xi,/22). 



(10) 



The Stillinger-David potential and its modified version 
are compared in Fig. l,b. The configuration of water 
molecules is depicted in Fig. l.a. 

3. Generalized Stillinger David Potential 

In this Section, the modified Stillinger-David potential 
is generalized further. The structure of the interaction 
potential between two water molecules (8) is adopted to 
remain invariant. Only the character of the interaction 
between oxygen atoms and hydrogen atoms in both the 
molecule and its neighbor changes considerably. In par- 
ticular, the form of the screening function 1 — i(r) is 
improved, and the fact is taken into account that the 
polarization of oxygen atoms depends, to some extent, 
on the electric field that emerges owing to a deformation 
of the electron shells of oxygen atoms induced by their 
immediate contact. 

3.1. Behavior of screening functions 

In the Stillinger-David approach, two screening func- 
tions, \ — K{r) and 1 — i(r), are used. We adopt that the 
function 1 — K{r) is equal to that obtained in work |12| : 



l-K{r)=r^/[r^ +F{r% (11) 
where 

F{r) = 1.855785223(r - Toh)^ exp[-8(r - ro^f] + 

-M6.95145727 exp[-2.702563425r] , 

and roH = 0.9584 A. 

The structure of the function \—L{r) is approximated, 
similarly to what was done in work |12| . by a combina- 
tion of an exponential function and a polynomial, 

l-L{r) = l-e-^"'-{l + Lir + L2r^+Lj.r^+Lir*). (12) 

However, the relevant coefficients are determined now, 
by using a new rigorous algorithm. It should be noted 
that nine parameters are to be determined simultane- 
ously; they include four coefficients Li, L2, L3, and L4 
in expression (12) and four parameters &i, pi, 62, and p2 
in formula (3). The listed parameters were determined 
under the following conditions: (i) the equilibrium dis- 
tance between the oxygen and hydrogen atoms in a water 



molecule was accepted to be Toh = 0.9584 A [4j; (ii) the 
angle between the directions from the oxygen atom to- 
ward hydrogen ones was taken 9 = 104.45° in the equilib- 
rium configuration (see work [4 ); (iii) the force constants 
= 2064.114 and #1^ = 91.5562 were taken from 
the molecular spectral characteristics (see work |19|); 
(iv) the final result had to reproduce the energy of inter- 
action between two water molecules that was obtained 
as a result of quantum-chemical calculations in work by 
Clementi and coauthors [TO^. We demanded that the 
function 1 — L{r) should change monotonously, because 
its nonmonotonicity would undesirably affect the behav- 
ior of the derivatives of the interaction energy of a water 
molecule. 

In addition, we adopted that the coefficient Lq — Li 
in the linear term r in the expansion of the function 
1 — L(r) in a power series of r in a vicinity of the point 
was equal to zero, i.e. Lq — Li. The coefficients Li, 
L2, and L4 were fitted to reproduce the value of force 
constant 

^ = ^^-^^^^ 

and to satisfy the condition 



do 



= 



(14) 



for the water molecule in the equilibrium configura- 
tion. Additionally, a requirement that the squared norm 
\\Li\\^ = Ll + L\ + Ll + Ll + L\ should be minimal 
was introduced into the determination procedure of the 
coefficients Li, L2, L3, and L4. A similar a priori con- 
dition was proposed by A.N. Tikhonov in work |20j to 
distinguish stable normal solutions of the systems of lin- 
ear equations. By definition, normal are those solutions, 
the moduli of which are close to zero. 

In such a way, we obtained that the coefficients Lq, 
Li, L2, ia, and L4 are 



Lq = 2.98, Li = 2.98, 



L2 = 0.92, L3 = 4.7044, L4 = 2.3580. 



(15) 



In what follows, the magnitudes of physical quantities 
are given to an accuracy of 5 significant figures, because 
it is to this accuracy that the experimental values of force 
constants were determined. 

The comparative behavior of the screening function 
1 — L{r) in the SD and GSD potentials and the corre- 
sponding function S{r) in the MPM potential is pre- 
sented in Fig. |2] We emphasize that the screening func- 
tion \ — L{r) is monotonous in the GSD potential (see 
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Fig. 2. Behavior of the screening functions in the SD, GSD, and 
MPM potentials 

formula (15)), in contrast to its behavior in the SD po- 
tential, and more adequately reproduces the extension 
of screening region in comparison with the function S{r) 
in the MPM potential g]. 

Note that the surface F{Lo, Li, L2, L^, L4) is very ir- 
regular in the multidimensional space of parameters. As 
a consequence, the equation ^. ^ — —23.133, which 
should have been used for the determination of the 
screening function coefficients Lqi -^ij -^2, -^3, and L4, 
turns out incompatible with Eqs. (13) and (14). Co- 
efficients (15) give the following value for this derivative: 
^ ^ — 205.9. For a similar reason, the experimental 

* was managed to be repro- 



dridr2 

value of the derivative 



dridr2 

duced in neither the SD potential nor the MPM one. 

For the determination of the parameters bi and pi, 
the following initial data were used: the aforesaid values 
(see Eq. (15)) for the coefficients Lq, Li, L2, L3, and L4 
of the screening function, the force constant 



92$ 



2 = 2064.114, 



(16) 



and the derivative 
dri 

Table 1. Parameters of the MPM, SD, and GSD 
approximations 



(17) 
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Fig. 3. Oxygen-hydrogen interaction energy in the SD, GSD, and 
MPM potentials 

which correspond to the equilibrium conffguration of a 
water molecule. The corresponding values obtained for 
the parameters of the GSD potential are 



hi = 3172.8, pi ^ 2.569. 



3.2. Parameters of the GSD Potential 



(18) 



In Table 1, the values for all parameters in the GSD 
potential are quoted, and a comparison with the corre- 
sponding values in the SD and MPM potentials is made. 

The values of force constants in the GSD potential, 
which were calculated for a water molecule with the help 
of the parameters presented in Table 1 are listed in Table 
2. For the sake of comparison. Table 2 also contains the 
values for the same constants in the SD and MPM poten- 
tials, as well as their values determined experimentally. 
One can see that the values of force constants obtained 
in the GSD potential are more exact than those in the 
SD potential. Concerning the MPM potential, the cor- 
responding force constants are almost identical to their 
experimental values. 

A considerable advantage of the GSD potential with 
respect to the SD and MPM ones becomes evident, when 
calculating the interaction energy between two water 
molecules (see below). 

A comparison of functions that describe the oxygen- 
hydrogen interaction within the MPM, SD, and GSD 

Table 2. Force constants for water molecule: experi- 
mental, MPM, SD, GSD 















ar2 




ariae 


dridr2 


Experimental |19j 


2064.114 


175.158 


91.556 


-23.133 


MPM 


2064.114 


175.158 


91.556 


286.94 


SD 


2064.114 


167.342 


-34.485 


117.57 


GSD 


2064.114 


167.342 


91.556 


205.9 



ISSN 2071-0186. Ukr. J. Phys. 2011. Vol. 56, No. 3 



229 



I.V. ZHYGANIUK 




Ron = 0.9584 A 



R 




dementi 
GSD ■ 



r, A 

Fig. 4. Relative arrangement of two water molecules used in work ^ig. 5. Comparison of the interaction energies between two water 

m 



molecules for the configuration depicted in Fig. 4 



potentials is made in Fig. 3. Figure 3 testifies that 
the core of an oxygen atom is too rigid in the MPM. 
This circumstance is very important, when calculating 
the interaction energy between two water molecules at 
distances that correspond to the dimer formation. 

3.3. General structure of the GSD potential 

The structure of the interaction potential between two 
water molecules in the GSD and MSD approximations 
is the same (see Eq. (8)). However, the calculation 
procedures for the dipole moment of oxygen in a water 
molecule and the parameters 62 and p2, which describe 
the influence of the neighbor water molecule on the oxy- 
gen dipole moment, change. The circumstance is taken 
into account that the electric field acting on the oxy- 
gen atom is a sum of the fields formed by the hydrogen 
atoms in the water molecule and the hydrogen and oxy- 
gen atoms belonging to the neighbor water molecules, on 
the one hand, and a component that emerges as a result 
of a deformation of the electron shells of oxygens, on the 
other hand. The latter effect is particularly important 
at distances between water molecules that correspond to 
the dimer formation. 

In accordance with all that, the dipole moment of an 
oxygen atom is determined by the formula 

^0lO2 

where the deformation field strength is 



j=3A 



Pi 



1 



^ e-P2»-Oi02 



O1O2 



P2 



0102 



(20) 



and Tr 



is the tensor of dipole- 

dipole interaction. The deformation field strength is a 
gradient of the repulsion force potential, E'''^-^ — — V$n. 

Hence, the dipole moments /Xo^ and /Xqj of oxygens 
are functions of the parameters (bi, pi) and (62, P2), 
respectively. The values of parameters bi and pi were 
determined above (see formula (18)). To find the pa- 
rameters 62 E^nd P2, let us calculate, using formula (8), 
the energy of interaction between two molecules in the 
configuration depicted in Fig. 4. 

For this configuration, five numerical values for the 
interaction energy at various distances between oxy- 
gen atoms were obtained in work [TU] on the basis of 
quantum-chemical calculations. The 62- and /32-values 
were so determined that the calculated curve optimally 
reproduced the positions of those points. The corre- 
sponding values obtained are 

62 = 42129.1, p2 = 2.59. 

A comparison between the energies calculated by for- 
mula (8) and reported in work jlOj is shown in Fig. 5. 
Below, when calculating the interaction energy for two 
water molecules, just those 62- and /92-values were used. 

In the MPM, the parameters 62 and p2 were deter- 
mined using the interaction energy <i>(o) — —8.47 be- 
tween water molecules in a dimer at the distance r^-^Q^ = 
2.96 A between oxygen atoms (see work |1]). 
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3.4- Dipole moment of an isolated water 
molecule 

The electric dipole moment of an isolated water molecule 
is determined as a sum of two antiparallel dipole-moment 
vectors, fi = /in + Mo- The dipole moment fi^ is 
defined by the spatial distribution of the centers of 
negative oxygen and positive hydrogen charges, /^h — 
gH(ri + ^2). The absolute value of dipole moment /Xh 
is = 2(7h''oh cos(|6') = 5.6281 D. The dipole moment 
/Xq of an oxygen atom emerges owing to the polarization 
of the electron shell of an oxygen anion by the electric 
fields of hydrogens in the water molecule. According to 
the results of work [12_, it equals 

Mo = -aq^ (^^ [1 ~ K{n)] + ^ [1 - ^('^2)]) . 

It is easy to calculate that /Iq — —3.7752 D. Therefore, 
the magnitude of the dipole moment fi is ^ = fi^ + fJ^o — 
1.8528 D. This value completely agrees with the absolute 
value of dipole moment of an isolated water molecule. 

3.5. Influence of a neighbor molecule on the 
dipole moment of a water molecule 

The modification of the dipole moment under the influ- 
ence of a neighbor molecule is one of the simplest man- 
ifestations of many-particle effects in the system. To 

estimate the influence of the second molecule, let us cal- 

(12) 

culate the ratio /io //^o between the dipole moment of 
an oxygen atom calculated in the pair approximation to 
the dipole moment of an isolated molecule. The depen- 
dence of the ratio /Iq //^o on the distance between the 
oxygen atoms of two neighbor molecules is exhibited in 
Fig. 6. 

4. Conclusions 

This work was devoted to the generalization of the 
known Stillinger-David potential <&sd, which is widely 
used for the description of the intermolecular interac- 
tion in water. Two versions were considered: the modi- 
fled Stillinger-David potential 'I'msd and the generalized 
Stillinger-David potential $gsd- The modifled poten- 
tial $MSD improves the behavior of the potential <i>sD 
at rather large distances between water molecules. The 
generalized potential "I'gsd takes additionally the polar- 
ization effects into account more adequately in compar- 
ison with the initial potential $sd- Moreover, the fact 
that the repulsion part of the interaction between oxy- 
gen atoms affects their polarization |21|[^ was taken 
into account. 




r. A 

f 12) 

Fig. 6. Dependence of the ratio fiQ //Jo on the distance between 
oxygen atoms in two water molecules 

The values obtained for the interaction energy be- 
tween water molecules with the use of the generalized 
Stillinger-David potential and in quantum-chemical cal- 
culations [10_ coincide within the whole range of inter- 
molecular distances. The results of detailed calculations 
of the ground state energy and the dimer vibrational fre- 
quencies on the basis of $msd and $gsd potentials will 
be presented in a separate work. 
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ySArAJIbHEHHH nOTEHIIIAJI CTIJlIHflXCEPA I flEBI^IA 
LB. }KusaHioK 
P 6 3 K) M e 

y poGoxi sanponoHOBaHO BflocKOHajieHHli nojiapH3au,iitHHit noTen- 
i;iaji CTijiiHfljKepa i /IsBifla fljia MijKMOjieKyjiapHoi BsaeMOflii' y bo- 
fli. C4)opMyjibOBaHo liTKHli ajiropHTM BHSHaHeHHa 4)yHKi(ii, aKa 
onHcye BsaeMOfliro OKcnreH— riflporen b MojieKyjii BOfln. Po3po5jie- 
Ho HOBHH niflxifl flo MOflejiKJBaHHa <|)yHKu,ii', m,o eKpanye sapafl- 
^iinojibHy B3ae]\io;j,iio iia MajiHx Bi^cxaHax. /Jjia npaBiijiBHoro onii- 
cy acHMnTOTHMHoi iiOBc;],iHKH MijKMOjieKyjiapHoro noTeHu,iajiy na 
flocTaTHbo BejiHKHx BiflCTaHax noTeHi;iaji CTijiiHfl»cepa i J\emf\a, 
3aBepmeHo B3aeMOflieio Mi>K ahhojibrhmh MOMSHTaMH oKCHremB. 
KpiM Toro, nojiapH3an;ifiHa CKjiaflOBa noTeHn;iajiy CTijiinflacepa i 
TleBifla flonoBHena flO^aHKOM, mo onHcye fle4>opMaii;iio ejieKTpoH- 
Hiix 060J10HOK OKCHreniB. y3arajiBHeHHa noTeHLi,iajiy CTijiin^jKepa 
i /leBi^a ^03B0jiae ycniiutio Bi^TBopHTH Bci ocHOBHi pe3yjibTaTH 
KBaHTOBO-xiMiHHHx po3paxyHKiB enepriii B3aeMOflii flBox MOjieKyji 

BOflH, OTpHMaHHX KjIBMeHTi. BHBHeHO nOBefllHKy flHnOJIBHOrO MO- 

MBHTy MOJieKyjiH BOflH aK 4'yHKi;ii MijKMOjieKyjiapHoT BiflCTani xa 
OTpHMaHO oii;iHKy HesBiflnnx flBOXHacTHHKOBHX e(J)eKTiB y BOfli. 
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